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Excitation of stellar pulsations 

Giinter Houdek 

Institute of Astronomy, University of Vienna, 1180 Vienna, Austria 

Abstract. In this review I present an overview of our current understanding of the 
physical mechanisms that are responsible for the excitation of pulsations in stars with 
surface convection zones. These are typically cooler stars such as the Delta Scuti stars, 
and stars supporting solar-like oscillations. 



1. Introduction 

Convection dynamics affects the driving and damping of pulsations particularly in stars 
with convectively unstable outer layers, including also the F stars which have rather 
shallow surface convection zones and in which solar-like oscillations have been ob- 
served. The type of stars in which the coupling between pulsation and the turbulent 
velocity field is important are, for example, the rapidly oscillating Ap stars, Delta Scuti 
stars, y Doradus, RR Lyrae stars, and stars which support stochastically excited oscil- 
lations. In this review I shall only address two of them, Delta Scuti stars and solar-type 
stars. 



2. Delta Scuti stars 

Delta Scuti (6 Set) stars have masses M between 1.5 < M / M Q < 2.5 and luminosities 
L between 0.6 < log L / L Q < 2.0. They are in the central or shell hydrogen burning 
phase with complex multiperiodic oscillation spectra, including both radial and nonra- 
dial modes (e.g., 79 pulsation modes have been detected in FG Vir, Breger et al. 2005). 
The oscillations are low-order p modes with relatively small amplitudes ranging from 
10~ 3 - O.lmmag and with periods between 18min and 8h. The modes are driven 
by the kappa mechanism in the second stage of helium ionization. Cooler 5 Set stars 
have substantial surface convection zones and, similarly to solar-like stars, the stability 
properties of the oscillation modes are crucially affected by the convection dynamics. In 
particular, the return to stability of low-order p modes at the cool boundary of the clas- 
sical instability strip is predicted only with the inclusion of convection dynamics in the 
stability analysis (e.g., Houdek etal. 1999, Houdek 2000, Xiong & Deng 2001, Dupret 
et al. 2005). These authors use different implementations for modelling the interaction 
of the turbulent velocity field with the pulsation. The currently most commonly used 
time-dependent convection formulations for studying the stability of stellar pulsations 
are those by Gough (1965, 1977a,b), Xiong (1977, 1989) and Grigahcene etal. (2005) 
the latter being a generalization of Unno's (1967) formulation. Although they all adopt 
the Boussinesq approximation to the fluid equations they differ substantially in detail. 
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A brief discussion of the main differences between these convection models was re- 
cently given by Houdek (2008). All three convection formulations model successfully 
the location of the red edge of the classical instability strip (IS), however, the physical 
mechanism responsible for the return to stability is very different in all three compu- 
tations. Dupret et al. (2005) report that it is predominantly the convective heat flux, 
Xiong & Deng (2001) the viscosity of the small-scale turbulence, and Houdek (2000) 
that it is predominantly the momentum flux (turbulent pressure p t ) that stabilizes the 
pulsation modes at the red edge of the IS. The results of these stability calculations 
are illustrated in Figure 1 in terms of accumulated work integrals W for stellar models 
located just outside the cool edge of the IS, i.e. for models for which the pulsation is 
found to be stable. 
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Figure 1. Accumulated work inte- 
grals of 5 Set models located just 
outside the red edge of the classical 
instability strip, i.e. the total work 
integral W <0. Panel (a) shows 
the results by Dupret et al. (2005) 
using the time-dependent mixing- 
length model based on Unno's 
(1967) and Gabriel's (1987, 1996) 
descriptions, and extended by Gri- 
gahceneetal.(2004). Panel (b) il- 
lustrates the results reported by 
Xiong & Deng (2007) using Xiong's 
(1989) Reynolds-stress-model-like time-dependent convection formulation. Panel (c) 
are the results by Houdek (2000) using Gough's (1977a,b) nonlocal, time-dependent 
mixing-length model. Panel (a) is adapted from Dupret et al. (2005), panel (b) from 
Xiong & Deng (2007), and (c) from Houdek (2000). 
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Nonadiabatic pulsation computations provide complex eigenfunctions and eigen- 
frequencies to = a> r + ia»; for modes with angular frequencies u r and growth/damping 
rates a)\. Work integrals W for complex eigenfrequencies were, for example, provided 
by Baker & Gough (1979) and can be divided into contributions arising from terms as- 
sociated with thermal energy equation (e.g., gas pressure perturbation, W g ), and terms 
associated with the momentum equation (e.g., momentum flux perturbation, W t ). 
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Figure 2. Linear stability coefficients ~ff := LoJu> r (positive values indicate in- 
stability) of the fundamental radial acoustic mode for an 1.7 M 6 Set star in the 
hydrogen-burning phase crossing the instability strip. Results are shown for two dif- 
ferent mixing-length parameters and for a calculation in which the effect of acoustic 
radiation in the equilibrium model was included according to expression (HJ. 



The upper left panel (a) of Figure 1 shows the results of Dupret et al. (2005). The 
individual contributions to the total W arising from the radiative flux, Wr, the convec- 
tive flux, W c , the turbulent pressure, W t (Wr + W c + W t ), and from the perturbation of 
the turbulent kinetic energy dissipation, W £ (Wr + W c + W £ ) are indicated by different 
line styles. It demonstrates the near cancellation effect between the contributions of the 
turbulent kinetic energy dissipation, W £ , and turbulent pressure, W t , making the contri- 
bution from the fluctuating convective heat flux, W c , the dominating damping term. 

The results of the calculations by Xiong & Deng (2001) are displayed in the upper 
right panel (b) of Figure 1. Contributions to the total W = W g + W t + W v arising 
from the gas pressure fluctuations, W g , the turbulent pressure fluctuations, W t , and from 
an effective viscosity from the small-scale turbulence, W v , are indicated by different 
line styles. The ionization zones of hydrogen (H) and helium (He) are indicated. In 
this calculation the dominating agent for making the pulsations stable is the damping 
contribution from the small-scale turbulence, W v . 

Panel (c) of Figure 1 shows the results of Houdek (2000) using Gough's (1977a,b) 
convection formulation. Contributions to W = W g + W t (solid curve) arising from the 
gas pressure perturbation, W g (dashed curve), and the turbulent pressure fluctuations, 
W t (dot-dashed curve), are indicated. The dotted curve is the ratio of the convective to 
the total heat flux, F c /F, and the ionization zones of H and He (5% to 95%) are indi- 
cated. In this calculation it is the contribution from the turbulent pressure fluctuation, 
W t , that stabilizes the pulsation mode. 
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It is interesting to note that all three convection descriptions include, although 
in different ways, the perturbations of the turbulent fluxes, but Gough (1977a), Xiong 
(1977, 1989) and Unnoetal. (1989) did not include the contribution W £ to the work 
integral because in the Boussinesq approximation (Spiegel & Veronis 1960) the viscous 
dissipation is neglected in the thermal energy equation. In practise, however, this term 
may be important. 

All three convection models need calibration for the mixing-length parameter a. 
The effect of varying a on the location of the IS for the fundamental radial acous- 
tic mode is illustrated in Figure |2] for an 1.7 M© star in the central hydrogen-burning 
phase. Decreasing a results in a significant shift of the cool boundary of the IS towards 
lower surface temperatures, whereas the blue edge is essentially unaffected. This comes 
about because reducing a diminishes the stabilizing effect of the turbulence on the pul- 
sations, which is, however, offset by the increasing efficacy of the surface convection 
with decreasing surface temperature T e g of the star. The observed location of the cool 
edge of the IS can therefore by used to calibrate a in 8 Set stars. 

Another effect, yet neglected in essentially all stability computations, is the ef- 
fect of acoustic radiation in the equilibrium model. The generation of acoustic waves 
is omitted in convection models which assume either the anelastic (Gough 1969) or 
Boussinesq approximation (Spiegel & Veronis 1960) to the fluid equations. However, 
through the generation of sound waves, kinetic energy from the turbulent motion will 
be converted into acoustic radiation (Lighthill 1952) and thus reduce the efficacy with 
which the motion might otherwise have released potential energy originating from the 
buoyancy forces. This effect may become important in stars with large turbulent Mach 
numbers M t := w/c, where w is the vertical rms component of the convective velocity 
field u = (u, v, w), and c is sound speed. In the phenomenological picture of an over- 
turning convective fluid element (eddy), which maintains balance between buoyancy 
forces and turbulent drag by continuous exchange of momentum with other elements 
and its surrounding (e.g., Unno 1967), the equation of motion for the turbulent element 
of vertical size i can be written as 

2w 2 _ , P ac 

= gaT' - -2 , (1) 
t l pw 

where T is the Eulerian temperature fluctuation, a the coefficient of thermal expansion, 
and g is the acceleration due to gravity. The rate of energy of acoustic radiation per unit 
volume, P ac is estimated according to the Lighthill-Proudman formula 

P ac=A ^!pr (2) 



€ 

where A and /u are the coefficients for emissivity and Mach-number dependence. If 
the emission of acoustic waves by homogeneous, isotropic turbulence is dominated by 
the largest convective eddies, the acoustic emission P ac scales with a Mach-number 
dependence /i = 5 (Lighthill 1952) and the emissivity coefficient A - 100 for a solar 
model (Stein 1968), which we also adopt here in our estimate for P ac . 

The effect of the acoustic radiation P ac in the mean model on the location of the 
instability strip for the fundamental radial acoustic mode is illustrated by the dot-dashed 
curve in Figure|2]for an 1.7 M Q star in the central hydrogen-burning phase. At the cool 
edge the effect of including P ac is similar to reducing the mixing length, however, at the 
hot (blue) edge of the IS P ac has a much larger effect than changing the vertical extent 
of the convective elements. 
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3. Solar-like stars 

Kepler has been providing data of solar-like oscillations with unprecedented quality in 
more than 2000 stars (e.g., Chaplin etal. 2011; Garcia etal. 2011; Huberetal. 2011b). 
One of the many important outcomes of the Kepler data analyses is the confirmation 
of the problem of reproducing theoretically the observed amplitudes of stochastically 
excited modes in stars hotter than the Sun. In general, Kepler observed photometric am- 
plitudes in stars located near the red edge of the IS that are rather similar to those found 
in the Sun, whereas theoretical amplitude estimates are up to 3-4 times larger than the 
Kepler values. This discrepancy in the oscillation amplitudes was already recognised 
before (Houdek 2006) when first ground-based spectroscopic observations were avail- 
able, such as for the F-star Procyon (e.g. Martic et al. 1999, Arentoft et al. 2008). It 
therefore behoves us to address the possible reasons why our current theory fails to 
reproduce observed amplitudes of solar-like oscillations in hotter stars. Reviews on the 
theory and modelling of stochastically excited modes were recently given, for example, 
by Appourchaux et al. (2009) and Houdek (2010). Therefore I shall summarize only the 
most important matters for our discussion. 

Solar-like oscillations are intrinsically damped but driven stochastically by the 
vigorous turbulent convection in the very outer stellar layers. The height H of a single 
peak in the observed oscillation power spectrum can be obtained from taking the Fourier 
transformation of a damped, harmonic oscillator for the surface displacement followed 
by an integration over frequency to obtain the total mean energy E in a single pulsation 
mode with (normalized) inertia / (e.g. Chaplin etal. 2005; Houdek 2006). Provided 
that the observing time is long compared to the mode lifetime the height H in units of 
cm 2 s -2 Hz -1 is given by 



where n is the damping rate or inverse of the mode lifetime, and P is the energy-supply 
rate in erg s -1 . If only the Reynolds stress driving term is considered in the equations of 
motions, the energy-supply rate P for radial models is given by (e.g. Chaplin et al. 2005) 



where R is surface radius, p t = (pww) is the (r, r) component of the (mean) Reynolds 
stress tensor (also known as turbulent pressure; angular brackets denote an ensemble 
average), % r is the normalized radial part of the pulsation eigenfunction and the 
product of cD and *F is a factor of unity accounting for the anisotropy of the turbu- 
lent velocity field u. The spectral function S accounts for contributions to P from the 
small-scale turbulence and includes the normalized spatial turbulent energy spectrum 
E(k) and the frequency-dependent factor Q(T£; <D r ); r\ :— A/ku^ is the correlation time 
scale of eddies with wave number k and velocity lit (A is a factor of order unity and ac- 
counts for uncertainties in defining r^), and k - kiln. For E(k) it has been common to 
adopt, for example, the Kolmogorov (Kolmogorov 1941) spectrum. For the frequency- 
dependent factor Q(Tjt;a>), however, which is used for evaluating the self-convolution 




(4) 



with 
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£l(u) T , Tk; r) = J Q.(co, Tk; r)£l(a> r -a>, r^; r) doj in equation ©, no satisfactory theory ex- 
ists. The two commonly adopted forms are 

- the Gaussian factor (Stein 1967) , 

n G (co;r k ) = ^ e -^^ 2 ; (6) 

V27T 

- the Lorentzian factor (Gough 1977b; Samadi et al. 2003; Chaplin et al. 2005) , 

n L (oj; T k ) = T ±_ 1 . - - (7) 

n V21n2 1 + (ojtJ V21n2) 2 



100.0 
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Figure 3. Maximum energy-supply rates P (equation!!]) of stochastically excited 
oscillations in 83 (symbols) solar-like stellar models, computed in the manner of 
Chaplin et al. (2005) using a Gaussian frequency factor (equation|6]l. The dashed line 
is a fit to the 83 model data by linear least squares, leading to an exponent s = 2.80 
in the scaling relation max(P)/max(P ) ~ (L/M) s , where L and M are in solar units 
(see also Table 1). 

The Lorentzian frequency factor is a result predicted for the largest, most-energetic 
eddies by the time-dependent mixing-length formulation of Gough (1977b), which 
decays more slowly with depth z and frequency a> than the Gaussian factor. Conse- 
quently a substantial fraction to the integrand of equation (0]) arises from eddies sit- 
uated in the deeper layers of the Sun, resulting in a larger acoustic excitation rate P. 
Samadi et al. (2003) reported that Stein & Nordlund's hydrodynamical simulations also 
suggest a Lorentzian frequency factor. However, Chaplin et al. (2005) reported that the 
Lorentzian time-correlation function leads to overestimated heights H at low frequen- 
cies for solar p modes. Recently, Belkacem et al. (2010) suggested to use a modified 
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Table 1 . Comparison of different stochastic excitation calculations. The correla- 
tion parameter A accounts for uncertainties in defining the correlation time scale of 
the turbulent eddies. The value of s is obtained from fitting {L/M) s (in solar units) 
to the computed values of max(P)/max(P ) of 83 solar-like stellar models by linear 
least squares (see Figure[3]). 
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Lorentzian frequency factor, by cutting off the high-frequency tail, i.e. setting it to 
zero, for frequencies o> > oje- The cutoff frequency u>e is associated with the inverse of 
the Eulerian microscale (Tennekes & Lumley 1972) and can be estimated with the help 
of Kaneda's (1993) 'random sweeping model approximation'. 

One of the many ongoing projects of the KASC (Kepler Asteroseismic Science 
Consortium) working group 1 "Solar-like p-mode oscillators" is the comparison of var- 
ious stochastic excitation models. In a preliminary exercise a grid of 83 solar-type 
main-sequence models was computed in the manner of Chaplin et al. (2005) from which 
energy-supply rates P were computed with the excitation models by Samadi et al. (2001, 
2003), using a Lorentzian frequency factor modified in the manner of Belkacem et al. 
(2010), and the excitation model by Chaplin et al. (2005) using a Gaussian frequency 
factor. All three computations assumed the same equilibrium models, pulsation eigen- 
functions and mode inertia. The maximum values of the energy-supply rate P for the 83 
solar-like stars with masses between 0.9 < M / M Q < 1.5, evaluated with the excitation 
model of Chaplin etal. (2005, see equation @]), are plotted in Figure |3] as a function of 
L/M, where L and M are in solar units. The maximum values of P scale with (L/M) s , 
where s = 2.80. The values for the exponent s for three different excitation model calcu- 
lations, all using the same equilibrium model and pulsation eigenfunctions are listed 
in Table 1. It interesting to note that the value of s varies only marginally between the 
different excitation models considered here, suggesting that the main problem of mod- 
elling the observed oscillation heights H (see expressions [3]-[5]) in hotter main-sequence 
stars is very likely related to a failure in modelling either the equilibrium structure (e.g. 
the convective velocity field and consequently p t ), the pulsation eigenfunctions or 
the mode damping rates n = -u)\. 

3.1. Pulsation eigenfunctions 

A very promising test of the pulsation theory, and in particular of the pulsation eigen- 
functions, independent of an excitation model, is provided by comparing estimated 
intensity-velocity amplitude ratios, AL/AV & , with observations (e.g. Houdek et al. 
1995, 1999; Jimenez et al. 1999; Jimenez 2002; Houdek 2009). The estimated ampli- 
tude ratios are obtained from taking the ratios of the complex pulsation eigenfunctions 
provided by the nonadiabatic stability computations, which must include the coupling 
of the turbulent velocity field with the oscillations. A comparison of solar observations 
(Schrijver et al. 1991) with model predictions are displayed in the left panel of Figure|4j 
where the model results are depicted for velocity amplitudes computed at different at- 
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Figure 4. Comparison of computed amplitude ratios of stochastically excited 
modes with observations. Left: results are shown for a solar model at different atmo- 
spheric heights (curves) and for the Sun (symbols). The thick solid curve is a running 
mean average of the observed solar data (adapted from Houdek et al. 1999). Right: 
results are shown for a model of Procyon A at different atmospheric heights (/i=300, 
600 km) and for two different atmospheres: a scaled VAL-C (Vernazzaet al. 1981) 
and an Eddington atmosphere. The solid curves (thin curves indicate 1 cr uncertain- 
ties) are the observations (Huber et al. 201 la). 



mospheric heights h. The square root of the mode kinetic energy per unit increment of 
radius r, which is proportional to rp 1 ^ 2 ^, increases rather slowly with height; the density 
p, however, decreases very rapidly and consequently the displacement eigenfunction g 
increases and the amplitude ratio decreases with height h. 

In the right panel of Figure|4]model results for the F5 star Procyon A are compared 
with observations (Guenther et al. 2008) from the MOST spacecraft and contempora- 
neous radial velocity measurements by Arentoft et al. (2008). Theoretical results are 
shown for two stellar atmospheres at two different atmospheric heights. For both stellar 
atmospheres the agreement with the observations is less satisfactory than in the solar 
case, indicating that we do not represent correctly the shape of the pulsation eigenfunc- 
tions in the outer stellar layers. Moreover, the modelled amplitude ratios significantly 
depend on the adopted atmospheric model. Consequently there is need for adopting 
more realistically computed atmospheres in the equilibrium models, particularly for 
stars with much higher surface temperatures than the Sun. 

Additionally to the effect of the atmospheric structure, the nonadiabatic eigenfunc- 
tions are also crucially modified by the dynamics of the turbulent velocity field in the 
convectively unstable surface layers. In particular the poorly modelled anisotropy of 
the turbulent velocity field u - (u, v, w) in the near-surface layers crucially affects the 
magnitude and frequency-dependence of the amplitude ratios (Houdek 201 1). In our 
current time-dependent convection model (Gough 1977a,b) the anisotropy O := u-u/w 2 
is parametrized by a constant value of order unity. Numerical simulations, however, 
indicate that O varies rapidly with height in the outer stellar layers, as illustrated in 
Figure [5] for models of the Sun (left panel) and Procyon A (right panel). It therefore 
behoves us, with the help of numerical simulations, to develop a model for the velocity 
anisotropy <I> for a more realistic description of the shapes of the convective cells as 
a function of stellar radius and consequently its effect on the pulsation eigenfunctions 
and frequencies. 
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Figure 5. Anisotropy O of the turbulent velocity field as a function of the depth 
variable pressure (/?) in the convectively unstable surface layers for models of the 
Sun (left panel) and Procyon A (right panel). Results are displayed for two different 
numerical simulations (red solid curve: R. Trampedach 1999, personal communi- 
cation; blue dashed curve: H.-G. Ludwig 2005, personal communication) and for 
our current version of a semi-analytical convection model of Gough (1977a,b; black 
horizontal line). Figure adopted from Houdek (201 1). 



3.2. Linear damping rates 

Another uncertainty in determining amplitudes of solar-like pulsations in hotter stars 
are the modelling of the mode lifetimes (Houdek 2006, Chaplin etal. 2009, Baudin et 
al. 201 1). Although an acceptable agreement between observed heights of solar low- 
degree modes in the Fourier spectrum and theoretical expectations has been achieved 
with the help of a time-dependent convection model (Houdek etal. 1999, Houdek & 
Gough 2002) the largest error in the theory lies in the calculation of the damping rates 
of low-order modes (Houdek etal. 2001, Chaplin etal. 2005, 2009). We located our 
principal problem in the stability analyses to be in the convection zone, below the upper 
boundary layer. An even larger discrepancy between observationally and theoretically 
inferred acoustic spectral linewidths was reported for solar-like stars that are hotter 
than the Sun (e.g., Houdek 2006, 2009). For example, with the help of hydrodynamical 
simulations (Stein et al. 2004), it was found that the theory underestimates the damping 
rates of the most prominent modes by a factor of two, and possibly more, for the F5 
solar-like star Procyon (Houdek 2006). This failure in the theory has recently been 
confirmed by using average linewidth measurements of solar-like oscillations in several 
hundreds of stars observed by Kepler (Chaplin 2010, personal communication, see also 
Appourchaux et al. 2012). 

Current theories (Chaplin et al. 2009) predict the average mode linewidth to scale 
with the fourth power of the stars effective temperature, whereas preliminary Kepler 
data suggest a steeper dependence, i.e. larger than four (see also Baudin etal. 2011 
using CoRoT, and Appourchaux et al. 2012 using Kepler data). As in the case of Pro- 
cyon, this indicates a failure of our current theory to predict mode linewidths in stars 
hotter than the Sun, most likely because of a physical mechanism still missing in our 
current theory. One such crucial mechanism is incoherent scattering at the inhomoge- 
neous upper boundary layer (Murray 1993; Goldreich & Murray 1994), which becomes 
increasingly more important for stars with higher masses and effective temperatures. 

Murray (1993) and Goldreich & Murray (1994) were one of the first to have 
stressed the importance of the contribution of incoherent scattering (77 sca tt) to the mode 
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Figure 6. Maximum values of the turbulent Mach number M t := w/c and tur- 
bulent pressure p l :- (pww) (angular brackets denote an ensemble average) versus 
stellar mass for models along the ZAMS (adapted from Houdek & Gough 1998). 



linewidths, and that it may dominate over other agents contributing to mode damping. 
They derived in the geometrical optics limit, at the top of the convective envelope, the 
scattering contribution 

a> r M t 2 

tyscatt : — — 77 (8) 

n(n + 1) 

for a radial mode with order n and frequency a> r , which is proportional to the squared 
turbulent Mach number M t . 

It was recognised by Houdek & Gough (1998) that the turbulent Mach number in- 
creases rapidly with stellar mass and surface temperature (see Figure ©. Consequently 
the scattering contribution to the mode linewidths may therefore dominate over the 
other contributions for stars with high surface temperatures, or more precisely, for stars 
with turbulent Mach numbers considerably larger than in the Sun. Hot stars have rela- 
tively shallow surface convection zones but very vigorous turbulent velocity fields (e.g. 
Houdek & Gough 1998). 

Goldreich & Murray (1994) adopted in their scattering model a time-independent, 
mixing-length-like, convection model and concluded that the mode energy is most ef- 
fectively scattered into other modes of similar frequency but higher spherical degree, 
with the surface gravity mode (f mode) being the ultimate recipient of the scattered 
energy. A time-dependent treatment of the turbulent velocity field will increase the fre- 
quency range of the modes over which the pulsation energy can be scattered (Gough 
1977a). An immediate question to be asked here is what the frequency distribution of 
this process will be, because the distribution determines not only the mode linewidths 
but also the shape of the spectral peaks in the oscillation power spectrum (e.g., Rast & 
Gough 1995; Rast 1999). Moreover, if the frequency behaviour of the damping pro- 
cess is similar to that of the stochastic excitation process, which is very likely the case 
(e.g. Houdek etal. 1999), we would be in the position to put further constraints on the 
frequency factor of the turbulent velocity spectrum in any stochastic excitation model 
(e.g. Houdek 2010). 
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